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We construct an analytical model for two channel, two-body scattering amplitudes, and then 
apply it in the description of the three-body J/tp — > K + K~tt° decay. In the construction of the 
partial wave amplitudes, we combine the low energy resonance region with the Regge asymptotic 
behavior determined from direct two-body production. We find that resonance production in the 
Kit channel in J/ip decays seems to differ from that observed in direct Kir production, while the 
mass distribution in the KK channel may be compatible. 
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I. INTRODUCTION 

Meson spectroscopy has played an important role 
in developing phenomenology and gaining insights into 
QCD in the non-perturbative domain. In an amplitude 
analysis of experimental data it is necessary to explore 
all of the available theoretical constraints because the ex- 
traction of resonance parameters requires the analysis of 
partial waves outside of the kinematic range of exper- 
imental data. In particular, amplitudes describing the 
mass distribution of a two-body subsystem in a quarko- 
nium decay may be different from those describing scat- 
tering of the same two particles. In this paper we focus 
on the isospin 1/2, 1 and spin-one, P-wave scattering am- 
plitudes in the 7T7t, KK and Kn channels and compare 
the phase shift data with two-body mass distributions 
from the three-body J/tp — > K + K~tt° decay. These am- 
plitudes are dominated by the ground state vector reso- 
nances p(770) and if* (892) that are well established as 
quark-antiquark, QCD bound states that are weakly cou- 
pled to the meson-meson continuum. There is also strong 
experimental evidence for higher mass vector resonances, 
although precisely how many and to what extent these 
are related to QCD single hadron states remains an open 
issue [THS]. 

In Table|T]we list the masses of the lowest vector meson 
states obtained from recent lattice QCD simulations [6 
and the quark potential model [7] and compare them to 
the data compiled by the Particle Data Group (PDG) [8]. 
Below 2 GeV the PDG lists two excited isovector reso- 
nances, the p'(1450) and p"(1700), that could have the 
quark model assignments of 25 and ID, respectively. In 
the lattice simulations of [6] the average pion mass is ap- 
proximately 400 MeV, which puts the p meson approxi- 
mately 130 MeV above its measured mass. Shifting the 
vector mesons masses from lattice computations down by 
130 MeV puts the first excited state around 1600 MeV, 
which is ~ 150 MeV higher than the measured mass of 
the p'(1450). While a resonance in the 1600 - 1700 MeV 
mass range can be clearly inferred from the tttt scatter- 
ing phase shift data [2], the experimental evidence for 
the p'(1450) is ambiguous [S]. The main motivation for 



the p'(1450) comes from the need to accommodate data 
on 47T production [5J[TU]. To the best of our knowledge, 
however, there has been no comprehensive analysis of all 
available P-wave data and the importance of the vari- 
ous inelastic channels, possibly even the dominant one, 
KK [5] is yet to be settled. For example, an alternative 
scenario that seems to be supported by the lattice results 
might be that the 2S and ID states are above 1.6 GeV 
while any residual strength corresponding to the PDG 
p'(1450) could due to residual interactions between pions 
and/or inelastic channel effects. 





P(l~ ) 


JT(l-) 


Lattice QCD [6] 


0.90 
1.8 


0.95 
1.8 


Quark Model [7] 


0.77(l 3 «Si) 
1.45(2 3 Si) 
1.66(l 3 L>i) 


0.90(l a Si) 
1.58(2 3 5i) 
1.78(l 3 £>i) 


PDG m 


0.775 
1.465 
1.720 


0.895 
1.414 
1.717 



TABLE I: Masses of the first few lowest-lying vector meson 
resonances. 



The vector mesons discussed above can be produced 
in J /if) — > KKn and 3n decays. In this paper we focus 
on the former; we studied the latter in Ref. ill . The 
J/tp — > KKir decay has been analyzed by the BESII 
Collaboration [12]. The Dalitz plot distribution of the 
K + K~ir° events has clearly visible sharp bands corre- 
sponding to the isospin-1/2, JsT ± *(892) and weaker bands 
in the first excited K'* resonance region. The distribu- 
tion is shown in Fig. [T] There is also a significant en- 
hancement in the low K + K~ invariant mass region. In 
the BESII analysis this broad band was associated with a 
new isovector P-wave resonance, the X(1570), decaying 
to K + K~ with the pole position at (1576 - 409i) MeV 
seen through a strong destructive interference with the 
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FIG. 1: The J/ip -> K + K~tt° Dalitz plot distribution from 
the BESII Collaboration [12]. 



ing [2J [TH1 E] ■ In [H] we successfully used these ampli- 
tudes to describe the J/ip —> ir + ir~ir° Dalitz distribution. 
In particular we have found that, since the //'(1700) is 
quite inelastic [2J, destructive interference with the vir- 
tual J/ip — > KKir — > 37T process is important in reducing 
the Dalitz plot intensity in the p" resonance region. We 
will investigate if it is possible that a similar phenomena 
is in operation in the KKn final state and whether the 
virtual J/ip — > 37r — > KKir decay may be responsible for 
the broad structure at low K + K~ invariant mass. 

This paper is organized as follows. The partial wave 
decomposition of the decay J/ip —> K + K~ir° is given in 
Sec. [ITJ In Sec. |III[ we discuss our P-wave amplitudes 
and compare with the BESII data of Fig. [I] Ideally, 
the set of partial waves that are developed here could be 
used in a full Dalitz plot analysis, but this requires a full 
knowledge of experimental acceptances and resolutions. 
In this work we simply compare, qualitatively, a sample of 
Dalitz plot distributions generated from our amplitudes 
with the BESII result of Fig. [TJ We include more details 
on the amplitude construction in the appendices. 



p(1700). There have been several theoretical attempts to 
explain this result [131 [T3] . 

In this work we address the following questions. Can 
the broad enhancement in the low-mass KK channel 
be described by the P-wave KK amplitudes determined 
from phase shift analysis? And, more generally, can the 
Dalitz plot distribution of KKir events in the J/ip decay 
be described in terms of KK and Kir amplitudes recon- 
structed from phase shift analysis? To do so, we use, 
and further develop (by incorporating asymptotic energy 
dependence), the P-wave 7T7r and KK amplitudes ini- 
tially constructed in [Til [15]. The amplitudes that we use 
have the correct analytical properties, satisfy two body 
unitarity and reproduce the known data on 7777 scatter- 



II. PARTIAL WAVE AMPLITUDES IN THE 

J/ip -> K+K'-k DECAY 

Denoting the four momenta by p±,o> P for K , ir° and 
J/ip respectively, the general expression for the J/ip — > 
K + K~ir° amplitude is given by, 



(tt°K + K- ,out\J/tp(\),in) = z(2tt) 4 5 4 ( E Pi - P T A 

J=o,± J 

(1) 

The Dalitz plot invariants are defined by = (pi +Pj) 2 
with i,j — ±,0 referring to K and the ir°, respectively. 
The general expression for the helicity amplitude of T\ 
is given by 



S,L M=±,0 



fj,,o( 



to)F+°(s +0 ) + Dl*(r-o)d 



)P Si u ( S -o)](2) 
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where N SL ^ = A /3(25' + 1)(S>; L0\lp) /4tt. Here A is the 
spin projection of the J/ip along the e + e~ beam axis, 
which together with x and y define a lab coordinate sys- 
tem, S is the spin of a two particle subsystem (the iso- 
bar), and L is the relative orbital angular momentum be- 
tween the isobar and the spectator meson. The rotation 



is given by three Euler angles, ry = 



''<./• 



which rotate the standard configuration in the (ij)k cou- 
pling scheme, to the actual one. In the standard config- 
uration of the (ij)k coupling J/ip is at rest, particle k 
has momentum along the negative z axis, and particles 



i and j have momenta in the xz plane with the parti- 
cle j moving in the positive x direction. The azimuthal 
and polar angles, (pij and "dij, are defined in the J/ip rest 
frame and refer to the actual direction of motion of the 
(ij) pair. Finally, ip\^ and 6\j are the azimuthal and the 
polar angle of the i-th particle in the {ij), two-particle 
(isobar) rest frame. 

The scalar form factors Fg L (sij) describe the dynamics 
of the decay in the isobar model i.e. under the assump- 
tion that in a given isobar channel the form factors are 
functions of the sub-energy of that isobar only. In the 
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L — S basis, the parity of the K + K~ir° state is given which we define as, 
by P = (— i) s '+ L + 1 and under charge conjugation the 

two isobar channels, \(K + tt°)K~) and \(K~ir°)K + ) are . s 

exchanged while the third isobar channel, \(K + K~)ir°) F+- = \ L-]pKK p+® — —Fa? = —F^f (3) 

is a charge conjugation eigenstate with the eigenvalue z 

(— l) s . Thus charge conjugation invariance implies that 

in Eq. ^ there are only two independent form factors and obtain, 



T a = E E N' SL ^[Dl*Jr + -)dl (6+_)F^(s + -)-D^ 

S,L ai=±1,0 

(4) 



with N' SL = Nslij,(1 + (— l) s+L )/2. The /i = com- further reduce the partial wave expansion to 
ponent vanishes due to parity conservation and we can 



Tx = ^N' SL1 {[Di^r + -) + Di% l (r + -)]dl (fft-) 

S,L 

- [UftCr+o) + D^ 1 (r +0 )]<o(^ )Fj f /( S+0 ) + [£>ft(r_„) + Di*_ 1 (r. )]df fi (9Z )F^(s. )}. (5) 

I 



Finally, it is useful to rewrite the above amplitude in 
terms of a single set of angles describing orientation of the 
decay plane. Using the relation between Euler rotations, 

r+_ = r_ r(0,x+,0)=r + or- 1 (0 ) 0,7r)r- 1 (0 1 x-,0) ) 

(6) 



where X+(X-) is the angle between K + (K~) and ir° and 
in the K + K~tt° rest frame enables to write T in terms 
of alone 



Tx =^N' SL1 [D% 1 (r + -) + Dl% 1 (r + -)][dl (et_) \^-F™{s + _) + df )O (0+ o )i^( S+o ) + df, (6r^F£f( s _ )]. 

S,L 

(7) 



The allowed quantum numbers in the K + K channel over the orientation of the decay plane and the other in- 

are S PC = 1 (p),3 (ps), and in the K ± tt° chan- dependent invariant, 
nels, S p = 1" (K* ) , 2+ {K* 2 ) , 3" {K* 3 ) , • • • . In the follow- 
ing we will assume that the Dalitz distribution can be sat- dT r»_o( s +-) 
urated with the lowest partial waves, i.e. P-wave in both j , = -^\A S H — / ds_o|T| , (8) 



II dy/S-\ 

K + K and K ± ir channels, and we test this hypothesis 
by studying the effect of the D-wave resonances in the Kir 
channels. Parity conservation implies S — L; therefore, 
in the following we will simply denote Fg L by F 1 ^ . The 1 — ( — l) s 

(unnormalized) J/ip partial decay width with respect to \T\ 2 = | N' SL1 [df (9+_) F^ K (s_| ) 

one of the Dalitz invariants (e.g. M K + K - = ^fs^ ) is ob- s,l 

tained by integrating the square of the decay amplitude _|_ ^+ ^p^ v ^ s+0 j -|_ d s )F^ 7V (s o)]| 2 (9) 



and 
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The overall normalization (N) is adjusted to match the 
measured number of events. It is \T\ 2 that determines 
the distribution of events in the Dalitz plot, (i.e. \T\ 2 = 
const, would give a flat distribution). The integration 



limits, s"g d "(s_| ) are roots of the equation which define 

the boundary of the Dalitz plot, 



s + s +0 S-o - (s +0 + 8-a){m%m 2 K + M 2 m 



2 k) s, 



. (m 



M 2 ml] 



•2(m 



K< 



M 2 m 



+ 2M 2 m 2 K ml) = 0. (10) 



Projections along M K +„ = ^/s +Q and M K -^ = ^/s_ 
axis can be defined analogously. 

In the following we discuss parameterizations of the 
form factors K and F^* in terms of two-body am- 
plitudes. Any parameters remaining in these parameter- 
izations, which are related to the production process as 
opposed to final state interactions should be determined 
by fitting the Dalitz distributions. As discussed in Sec. [I] 
we do not fit the published Dalitz distribution, but in- 
stead show the predicted distributions for specific values 
of these parameters. 



III. THEORETICAL MODEL FOR FORM 
FACTORS 



Unitarity relates production form factors to two-body 
amplitudes. In [TTJ [15] we constructed analytical repre- 
sentations for the isovector, P- wave, two-body, nir and 
KK amplitudes. Here we further extend the analysis of 
[llj by constraining the high energy behavior, and ex- 
tend the approach to the Kir channel. We begin with 
a if -matrix, phenomenological parameterization of the 
known data (on the real axis) on phase shifts and elas- 
ticities. Even though the if -matrix offers an analytical 
representation for the amplitude, it often leads to spuri- 
ous poles and zeros of the amplitude when extrapolated 
outside the physical region. Therefore we use the analyti- 
cal representation for phase shifts and inelasticity via the 
if -matrix only in the data region and smoothly extrapo- 
late to match with the asymptotic behavior of the partial 
waves at high energies. We then use the amplitudes con- 
structed this way over the whole physical energy range as 
input into the Omnes-Muskhelishvili integral to construct 
the part of the scattering amplitude regular on the left 
side of the complex s-plane. With the N(s)/D(s) rep- 
resentation, which is described below, we determine the 
amplitude over the entire s-plane. Finally we solve the 
unitarity relation for the form factors and write the J/tp 
decay amplitude in terms of the denominator functions 
D(s) and production functions c Q (s). In the following we 
describe these steps in a little more detail. All details of 
the amplitude construction are given in the Appendix. 



A. Amplitude Parameterization 

In [TT] to describe the high energy limit of the isovector 
P-wave, the following hypothesis was made: the S matrix 
is saturated by two channels, tttt and if if and the elastic 
channel phase shifts asymptotically approach a multiple 
of 7r with elasticity rj approaching 1. Even though J/tp 
decays probe only a limited energy range, and are quite 
insensitive to details of the asymptotic behavior we might 
as well use a different hypothesis that is better rooted in 
high energy phenomenology. It is known that at high 
energies, elastic cross sections slowly grow with energy 
almost approaching the Froissart bound. This implies 
that at impact parameter larger than the interaction re- 
gion 0(1 fm) there is no interaction while the low par- 
tial waves are suppressed as if scattering from a "gray 
disk." The low partial waves correspond to L << Lq(s) 
where L (s) ~ Vs/2 fm and while the interaction radius 
grows logarithmically with energy, the scattering of the 
low partial waves becomes logarithmically suppressed, 
i.e. t]l ~ 1 — 0(1/ logs). In the language of Regge 
exchanges this picture corresponds to the Pomeron ex- 
change at high energies. Furthermore, since asymptoti- 
cally the number of inelastic channels grows rapidly, each 
individual inelastic amplitude, e.g itit — > KK is expected 
to fall off with energy, and is represented by exchange of 
non- vacuum quantum numbers, aka meson Regge trajec- 
tories. The hypothesis of two-channel dominance in the 
high energy limit is therefore not necessarily well justified 
and in the following we adopt the Regge picture of high 
energy scattering. Matching the if -matrix parameteri- 
zation of the low energy data with Regge asymptotics, 
leads to amplitudes of the form (we drop the angular 
momentum label on the partial wave), 



4-Krnatrix ( \ „ ^ „ 
aJ3 y ' ' °low 

iHeqqe / \ . 

*«/ ( s )> s > s high 



(11) 



with tK™ atnx (s) and t* e p " e {s) determined from K- 
matrix fits to the low energy data and Regge fits to the 
high energy fixed £-data, respectively. Greek indices de- 
note two body channels, i.e. a — (i,j) — tttt, KK, etc. 
For energies between si ow and Shigh, we smoothly connect 
both real and imaginary parts of the if -matrix and Regge 
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amplitudes. The denominator function in the N/D pa- 
rameterization 



t a p(s) 



(12) 



is then obtained from the phase of the scattering am- 
plitude using the Omnes-Muskhelishvili solution of the 
unitarity relation (s t h > min(s a , sp) where s a is the a- 
channel threshold) 



lmD a p(s) 

and is given by 

D a p{s) 



-i f ds'- 



(13) 



(14) 



where we conveniently normalized D a p(0) — 1. The nu- 
merator functions N a p(s) are given by the largely un- 
known discontinuity of the amplitudes on the left hand 
cut. For the purpose of solving the unitarity relation 
for the J/ip decay form factors, which will be discussed 
below (c/. Eq. (16l), it is convenient to have iV a ^(s)'s 



for all intervening a, (3 channels having the same analyt- 
ical form. This is certainly a simplifying approximation, 
nevertheless we have found that with a simple parame- 
terization 



A 



a0 



SL 



(15) 



and 
Eqs. 



with the two-body amplitudes given by 
(12 1,(14), it is indeed possible to obtain good 

. phase shifts 



fits to the two-body scattering data, i 
and elasticity. 

Having constructed the two-body amplitudes, the next 
step is to relate them to the production form factors. 
This is done through the unitarity relations which relate 
the imaginary part of the form factors to the two-body 
amplitudes, 







(16) 



with a pt representing the elastic L-partial wave scatter- 
ing amplitude between two-body channels a = (ij) and 
(3 — (i'j 1 )- F is the reduced form factor (with the barrier 
factor removed), 

FH( S )= q ? j (s)p%( S )FV(s), (17) 

with Qij(s) being the relative momentum between mesons 
i and j, 



qij(s) 



[s - (m, + TOj) 2 ][s - (m l 



As 



and 



Pk{s) 



s- (M + m k ) 2 ][s- (M-m k ) 2 ] 



AM 2 



(18) 



(19) 



with M being the J/ip mass, the relative momentum be- 
tween the (ij) pair and the spectator meson k. p a (s) — 
^Qij/Vs describes the two-particle phase space. It is 
straightforward to show that if the scattering amplitude 
is dominated by a single resonance below inelastic thresh- 
old (p — p a , — for j3 ^ a) the solution of the uni- 
tarity condition for F is 



F£(s)=c(s)BW^( S ) 



c(s) 

- im R T L (s)' 



(20) 



where BW^(s) is the Breit-Wigner amplitude (with an 
energy dependent width T^(s)) and c(s) is a real polyno- 
mial in s. In the general multiple-channel case, with the 
two-body amplitudes all parameterized by the same nu- 
merator function, as in Eq. ( 15 ) the solution to Eq. ( 16 ) 
is given by [T5] . 



F a (s) 



\ - cp(s) 

2^ Daf}f 8 ) 



(21) 



with cg(s) being analytic functions in the right hand 
plane and Imc Q (s) = for s > 0. 



B. Results 

As discussed in Sec. |T] the original BESII analysis was 
based on the isobar, resonance parameterization of all 
three two-body channels. In the absence of a known 
isovector P-wave KK resonance to describe the low mass 
KK enhancement, it was necessary to introduce a new 
resonance, the X(1570). The tttt phase isovector P- 
wave shift data, however, points to significant inelasticity 
above 1.6 GeV, which following [2] we have attributed to 
the KK channel. The effect of the coupled irw and KK 
channels on the K + K~ 



mass distribution which fol- 
lows from Eq. (21 1 is shown in Fig. [2j 

In Fig. [3] we show the Dalitz distribution obtained us- 
ing the single Kir channel amplitude (details discussed 
in the Appendix). Besides the if* (892) peaks, bands at 
Mjin = 1-75 GeV are clearly visible in both K + tt° and 
if~7r° mass projections. These are due to the if* (1680) 
resonance clearly seen in the Kir phase shift analysis [T§1 - 
I21j but apparently not so in the Kix production from 
the J/ip decay (c/. Fig. [T]). This clear discrepancy in- 
dicates that it is not sufficient to use a single channel 
Ktt amplitude in the parameterization of the correspond- 
ing form factor in the J /ip decay. As discussed in the 
Appendix the Kir amplitude is indeed inelastic above 
Mktt ~ 1.5 GeV with a possibility of a large coupling to 
the K*(892)n channel. 

Finally, in Fig. [2] we show the Dalitz distribution ob- 
tained with a combination of three amplitudes, KK — > 
KK, tttt — » KK and Ktt — > Ktt with relative production 
coefficients, c a (s) chosen to best match the observed dis- 
tribution in Fig. |l]). While the low mass KK region 
seems to be fairly well described, the resonance struc- 
tures in he Ktt channel do not match between the elastic 
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FIG. 2: Dalitz plot distribution obtained using, in Eq. ( |21[ ), 
a single two-body KK — > KK amplitude (top) and a single 
two-body 7T7T — > KK amplitude (bottom) (i.e. with c K ^ = 1 
(ctttt = 1) for the top (bottom) and c a (s) = for all other 
waves) . 




FIG. 3: As in Fig. [2] but with a single two-body, Kty — > Kir 
amplitude. 




FIG. 4: As in Fig. [5] but with three amplitudes, KK 
KK, im — > KK and -KV — > Kn with relative production 
coefficients satisfying cktt{s) : c Jf ^-(s) : c^^(s) = 1 : 0.3 : 
-0.7. 



tKn-^Kir and J/V> decay amplitude. The irir — > if AT and 
— > KK amplitudes behave rather smoothly in the 
region corresponding to the Kir resonances and do not 
give enough strength to reducing the peak from the sec- 
ond K* resonance region. Thus we anticipate that the 
discrepancy is due to inelasticities in the Kir channel it- 
self. Since we are only comparing Dalitz distributions 
as opposed to fitting data, we do not attempt to fur- 
ther improve the comparison. It is worth noting that the 
fsT*(1410) listed in the PDG is indeed quite inelastic with 
only a 6.6% branching to Kir. 



IV. DISCUSSION AND CONCLUSION 

Based on unitarity and analyticity we have constructed 
a set of analytical two-body amplitudes, which imple- 
ment the known phase shift data. These extend our 
previous work in coupled channel P-wave 7T7t and KK 
systems and the J/ip — > 3ir decay [11]. The two-body 
amplitudes are only an approximation to the three-body 
decay, nevertheless they provide a useful starting point 
and should match below inelastic thresholds. We com- 
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FIG. 5: Phase shift of the P-wave tttt amplitude. Data is 
taken from [2] (circles) ,[16] (triangles) , and |17| (squares). 
The solid line is the result of the fit to 5^ and r\ with the 
analytical if-matrix representation described in the text. 
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FIG. 6: Same as in Fig. [5] for the inelasticity 77 



Appendix A: Analytical model for the P-wave 
isovector nn — > nn, nn — s> KK and KK — > KK 
amplitudes 



pared the analysis of the J/ip decay with these ampli- 
tudes to the original analysis of the BESII collaboration, 
which was based on the isobar model with coherent Brcit- 
Wigner resonances. The isobar model with the known, 
low mass resonances only and without inelasticities can- 
not faithfully produce the broad structure of low K + K~ 
invariant mass, which is why in the BESII analysis an 
additional P-wave resonance A r (1576) coupled to K + K~ 
was introduced. Our preliminary study indicates that 
the KK low-mass region may be described by the inelas- 
ticity in the irir — > 7T7T wave if attributed to the coupling 
between irir and KK channels. A single Kit — > Kir am- 
plitude is strongly affected by the second vector K* (1680) 
resonance as observed in Kn phase shift analysis. How- 
ever, in J/tp decay this resonance seems to be suppressed. 
It is worth noting that a similar suppression of the first 
excited isovector-vector resonance is also observed in the 
3n decay of J/tp [IT] . 
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1. if-matrix parameterization, (s < si ow ) 

We use a two channel if -matrix .11] to fit the data on 
7T7T — > 7T7T P -wav e phase shift and elasticity 77 from [2J 
Hl)l IT7] (Fig. 5j6 1 . With the S'-matrix saturated by two 
channels, the model makes a prediction for the phase 
shift in the KK — »• KK channel. In this section a, (3 = 
7r, K correspond to the two body channels 7T7t and KK, 
respectively. The 2-channcl if-matrix representation is 
given by 

[* _1 0)W = [K-^s)]^ + S aP (s - s a )I a {s), (Al) 



where 

I a (s) = I a (0) 



ds' \ / 1 



(A2) 



A convenient choice for the subtraction constant, / Q (0), 
is to take KeI a (Mp) — so that one of the poles of K V7r 
corresponds to the Breit-Wigner mass squared, = 
(0.77 GeV) 2 , of the p meson. Using the general two-pole 
parameterization of the K matrix, 



K = '■ 

™ Ml 



KirK — Kktt — 



PI 



s s 2 - 

ft IT fin 

s 2 - s 



+ 7tt 



TttK) 



K K k = 



S2 



IKK 



(A3) 



where a 2 . = r p M 2 /(M 2 - s T ) 3 / 2 and fitting the P-wave 
7T7T phase shift, 5 n , and the elasticity, rj, we obtain T p = 
0.140 GeV, and 



/s 2 = 1.4708 GeV, 
7™ = 5.62 x 10~ 2 , 



Pit = 0.199, p K = 0.899, 
lnK = 0.104, lKK = 1.525, 

(A4) 
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with the 7's in units of GeV - . The comparison of the 
phase shift and the inelasticity obtained with this pa- 
rameterization with the data is shown in Fig. |5|6| To 
illustrate unphysical features of the if -matrix parame- 



terization we rewrite Eq. ( Al I using the standard N/D 
representation for t a p = t a p/(4q a qp). With the normal- 
ization D a p(0) = 1 we obtain, 



N^(s) = A 
N K k(s) = \ KK 



r, D n7r (s) = exp / ds -— r , 

{s - s LA ){s - s L>2 ) V n Js„ s'{s'~s)J 

r, D nK (s) = ^ r exp / 

n ( \ ( S f J ' ^KK(s') \ 



s - z KK 
{s - s L ,i)(s - s L . 2 ) 



ds' 



s'{s' - s)J ' 



(A5) 



with A, 



5.649, \ KK = 2.271 and K K = 3.048 GeV 



In this if-matrix model, the left hand cut of N is re- 



duced to two poles at sl, 



-13.87GeV^ and s L 



-0.787 GeV , respectively. There are also first or- 



der zeros in N a p at z. 



-0.867 GeV^ and z KK 



-13.78GeV 2 . Above the KK threshold the phase of the 
inelastic amplitude (p^K is given by 4> v k = ^it+^k- From 
the K matrix we find that, asymptotically, (jj-xK (00) — 
2ir, which corresponds to two CDD poles: one at the 
p mass, si.Trif = Mp, and the other at S2,-kK = ^2 + 
PkPk/i-kK — 3.884 GeV 2 . Thus, while the K matrix pa- 
rameterization faithfully reproduces the rnr phase shift 
and elasticity in the whole available energy range, from 
7T7T threshold up to 1.9 GeV, extrapolation beyond this 
range is problematic. The rapid decrease of (j)^^ around 
s ~ 6 GeV 2 seems unphysical. In the ttt: — > KK channel, 
the two CDD poles at m 2 and s 2 + PwPk/JttK are clearly 
an artifact of the pole parameterization of the iVmatrix. 
A CDD pole in the inelastic channel above threshold (c/. 
the pole at S2^k = 3.884 GeV 2 ) leads to a discontinuity 
in a phase shift and is unphysical. It also implies van- 
ishing inelasticity, rj = 1 at this energy. A pole between 
7T7T and KK thresholds is admissible, e.g. the pole at 
sx,irK = w 2 , but its strict overlap with the p mass is 
also an artifact of the parameterization. Since the phase 
space available in J/ip decay extends up to s T7r ~ 9GeV 2 
we need to remove these unphysical features of the K- 
matrix amplitude. As discussed in Sec. [ill] we do this 
by using the iVmatrix amplitudes below si ow and above 
Shigh we will use Regge parameterization. 

2. Regge parameterization (s > s^igh) 

Regge analysis of irir — > irn scattering has been stud- 
ied recently in [22H24] and here we use the results of 
|24j . Parameters in Regge amplitudes were constrained 
by analyzing NN, irN and tttt scattering data. For com- 
pleteness we give the relevant formulas below. 

• 7T7T — > 7T7T 



Regge parameterization involves the Regge poles corre- 
sponding to i-channel exchange of the Pomeron(P), the 
P' (associated with the /2(1270) trajectory) and the p. 
The i-channcl isospin amplitudes are given by 

F¥r o \t,s,u) 

1 1 e -iira P {t) 1 1 _— i-K<x pl {t) 
= — 77^P(s,t)-— 7}S -P , (8,t), 

smirap(t) smirap'(t) 

(A6) 



1 _ e -ma p (t) l + ap (t) 

sin7ra p (i) p l + a p (0) 



[\ + d p t]e bt {s/s) a ^\ 



ImP^= 2 )(t, s, u) = /3 2 e bt (s/s) a ^ t)+a ^- 1 . 



where (s = 1 GeV) 
P(s,t) = (3 P a P (t) 



1 + ap ® e*( a /sr*n, 



(A7) 



(A8) 



(A9) 



P'(s,t) = 



ap,(t)[l + a P ,(t)} 



'(^""^.(AlO) 



a P ,(0)[l + a P ,(0)} 
and the trajectories are given by 
ap(t) = ap(0) + ta'p, 

a P ,{t) = a p {t) = a p (0) + ta' p + ^t 2 a" p . (All) 

Numerical values of all parameters are given in Eqs. (B5), 
(B6) of [23] • The s-channel isospin, partial wave ampli- 
tudes are normalized according to 

*S } M^) = (y2)' T ^(2L + l)^ ) ( S )P i (cos0), 

L 

(A12) 

where (V2) a is the identical particle symmetry factor: 
cr = 2 for 7T7T <H- 7T7T, (T = 1 for 7r7r <H- KK and cr = for 
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KK <H- KK. The s-channel amplitudes with I s = 0, 2 
are symmetric under t ■<+ u exchange, and the I s = 1 
amplitude is antisymmetric and s +>• t crossing leads to 
the following relation between the s and the i-channcl 
isospin amplitudes, 



(> 



F^(t,s,u)-(t^u). (A13) 



The (t O u) exchange brings in the u-channel Regge 
poles (these were ignored in [21] where only the forward 
t = limit was considered). Finally, projecting out the 
P-wave amplitude yields, 



(dcos9) cos( 



x [±F£=°\t,s,u) + ±Fl I r 1 \t,s,u) 

- ^= 2 )(t, a> tt)-(t-^u)]. (A14) 

The angular integration is done numerically. The leading 
asymptotic behavior due to Pomeron exchange can be 
calculated analytically and is given by, 

e 99e (*)^~* j\dcos9)cos9[P( S ,t)-P( S ,u)} 

In ,Q \ 

(A15) 



. 7T -3a P + 2(6 + a p lns) 

l ~ ,PP /l , _/ i_ \ 2 s 



24 



(6 + a'p In s) 



To combine the If -matrix (s < s; ou) ) with the Regge 
projected (s > Shigh) amplitudes into the full P-wave 
7T7T — > irir amplitude, 



l-Kmatrix 



(s),s < s low 
t^{s),s>s high 



(A16) 



^77(5) - 

we choose ^/S[ ow — 2.20 GeV and ^/Shigh = 2.56 GeV, 
and use a simple analytical formula to smoothly join the 
two amplitudes between si ow and s^igh- The result is 
shown in Fig. [7j 

• 7T7T — > if_K^ 

Asymptotically the i-channel amplitude it is dominated 
by the K* trajectory, 



F 



K 



(t,s,u) 



suiTro^) 2a K *(0) + l 



(A17) 



Following 



we use & = 2.4 Get/ 2 , and OK-«(i) 



0.352 + a' K »i with a^. = 0.882GeF- 2 . The s-channel 
I s = 1/2 amplitude is antisymmetric under t o u ex- 
change, 



p 



K 1 



s, f, u) 



-F 



K 



(s, u, t). 



(A18) 





FIG. 7: Real (top) and imaginary (bottom) parts of the 
isovector, P- weave amplitude, t 7T7T (s) (solid lines). Dashed 
line is the result of the /^-matrix parameterization. 



and from sf>i crossing we obtain, 



(s,t,u) 



(A19) 

In terms of P^ (S; i, u) the properly normalized P- 
wave in ttt: — > KK is finally given by 



(dcos 0) cos 9 



3 Kk 



(t, s, u) — (i — > u) 



(A20) 



We can fix /3k* by matching our formula in Eq. (A17) 



to Eq. (81) in [23J in the limit of t — > (forward direc 
tion). Taking into account differences in normalization 
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employed here and used in we find 



ImF„ 



ttK 



(t,S,u) I g-K»,t-K) 

3 A 



with A = 1.82 taken from [25] and 
3 A 



P. 



K* 



47rr[ajf.(0)] 



(a' K ,s) 



a K ,(0) 



(A21) 



4^T[a^(0)] 
Asymptotically, t^^ 9e (s) approaches 



= 0.172. 



(A22) 



t 



Regge 
K 



I _ e -iira K * (0) 

sin ItOLK* (0) 
I _ e -«7rajf* (0) 



^8WJ dcosg)cos ^(o) + r 



-> u) 



(3 K *a' K , (1 + 2a*. (0)) [6 + a' K , ln(a' K ,s)} - 2a' K , 



sin to x .(0) 3^ 2a K * (0) + 1 



[6 + aJf. ln(a^.s)F 



(A23) 



]The complete amplitude is given by, 



tKmatrix 
.Regge 



(s), S < S tow 
t Jxegge^^ s > g ^ 



(A24) 



where we choose ^/si ow = 2.5GeV and ^/Shigh — 3GeV 
and it is shown in Fig. [8j 

• KK -> KK 

Asymptotically we only retain the Pomeron exchange, 



F 



KK 



(t,s,u) 



1 -I- p-ina.p{t) 

— ■ T^rPp a P\ t ) 

sin na p (t) 



l + a P (t) 



e bt {s/s) apW . 

(A25) 

with atp(t) = ap(0) + ta' P and all other parameters, ex- 
cept Pp K taken from [22] [24], while for the Pomeron 

coupling to KK we use relation = { f -f^) 2 {f ( / ] ) 2 = 

1.15, where the values of -jfpj and jip — (/i ) are taken 
from [22] [24]. From sttt crossing, 



F%p\s,t,u) = \F^\t,s,u) - l -F^\ t ,s,u), 



(A26) 



For the Pomeron contribution to the s-channel P-wave 
we thus find 

el 



Regge 
KK 



7T 



,1 



(dcos0) cos6^F^ (t, s,u). (A27) 



-l 



Asymptotically, t I ^ e ^ 9e (s) is given by 



t 



Regge 
KK 

n 

t aKK-- 3 ^ +2(6 + a' P lns) nap(0 )_i 



/ X (d cos 0) cos gap (t) 1 + " P W e bt W 



(6 + a' p In s) 2 



(A28) 



In the full amplitude, 



6 if if 



4-Kmatrix 



(5), 5 < S k 



^KK (^)' ^ ' > ^high 



(A29) 



we take y/si ow = 1-62 GeV and JshAgh = 3 GcV for 
real parts of amplitudes, and ^fsi ow — 1.64 GeV and 
y/Shigh = 1-8 GcV for imaginary parts of amplitudes. 
The different choice for the real and imaginary parts al- 
lows for a smoother connection with the Regge asymp- 
totics. The phase of t^f9. 9e (s) asymptotically approaches 
7r/2 but the phase of t^R atnx (s) has a sharp drop above 



1.65 GeV (see right plot in Fig. 10 1. Therefore, choosing 
^Jsi ow ~ 1.64 GeV allows for a continuous match between 
the phases of t K g(s), as show in Fig. [9] 



3. Phases of amplitudes and D functions 

From Regge parameterizations we find the following 
asymptotic behavior for the phases <f> a f}{ s ) of the com- 
plete amplitudes, 



arctan 



sin 7rap(0) 



1 + cos7rap(0) 2 ' 

, . siniraK* (0) , 
birK -^2tt + arctan — - : 

1 — COS -KOIK' (0) 



2tt 



b K j{ —> arctan[— 



sin7rap(0) . 
1 + cos 7rap(0) 



7T 

2' 



(A30) 



r, (A31) 
(A32) 



These are shown in Fig. 10 For the D functions test lead 
to the following asymptotic limits (c/. Eq. (14)) 



1 



1 



D nw (s) s i ' D 7rK (s) 
1 i 



d kk(s) S5 



I _ e -iira K , (0) -y 

sin7ra^* (0) s 2+ s ' 

(A33) 
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1 



5 10 



FIG. 8: Real (top) and imaginary (bottom) parts of the 
isovector, P-wave amplitude, t n K(s)/ (qvgif) (solid lines). 
The dashed line is the result of the A'-matrix parameteri- 
zation. 




■ (GeV ! ) 



FIG. 9: Real (top) and imaginary (bottom) parts of the 
isovector, P-wave amplitude, t K R ( s ) (solid lines). The dashed 
line is the result of the if-matrix parameterization. 



and 



Appendix B: Analytical model for the P-wave 
ttK — > irK amplitude 

1. if-matrix parameterization, (s < si ow ) 

To fit the phase shift data on ttK scattering we use 
a two-channel if-matrix model, with the two channels 
being Ktt and K* (892)-7r, and in the second channel treat 
the K* as a stable particle, (i.e. we ignore cuts on the 
third sheet). Similarly to the 7T7t, KK case for the K- 
matrix representation of Ktt and K*(892)n amplitudes 
we write 

(Bl) 



where i afj = t a p/{Aq a qp), 



is Sq)(s Sq) ± o 

q a = \j — , s a = [m a ± , 



4s 

mi = m K , m 2 = M K , 



(892): 



(B2) 



I a (s) = I a (0) 



, l( S , 3 ■ (B3) 

7T S (S — S) 



A convenient choice for the subtraction constant, i Q (0), 
is to take Rei Q (M^.» , g92 s ) = so that one of the poles 

of if ii is located at mass squared of the if* (892), m\. 
In terms of phase shift and inelasticity the Ktt and 
if*(892)7r amplitudes are given by 



til — i*22 — 



*12 — *21 



2ipi " 2ip 2 

y/l — ^2 e i(<5n+<5 2 2) 

2VPiP2 



(B4) 



where p a {s) = \J (1 — s f-)(l — ^f-)- The denominator 
D a/ 3 of the Ktt and if*(892)7r amplitudes are defined by 
Omnes-Muskhclishvili function 



D a p(s) = cxp -- / 



ds'4^). (B5) 
2 s'(s' — S) I 



12 





FIG. 10: Phase of the tttt (top), tyK (middle) and KK 
(bottom) P-wave amplitude. The dashed line is the result of 
the if-matrix parameterization from Eq. |A5|. 



To fit the P-wave phase shift data [TT?H2"T] we use a three- 
pole parameterization of the K-matrix 



K n = 



K22 

K 12 
where 



M 2 



s 2 - s 



S3 



, (0) 
+ 7ii 



\2 



S 2 - S 
K 2 1 



S3 ~ 



S2 



, (0) , (1) 
+ 722 + 722 S. 

AlA 2 (0) (1) 

+7l2 + 7l2 S : 

S3 — s 



V (892)^|:»( 8 92) 



[(M 2 



if* (892) 



*(892) 



7^, 



(B6) 



(B7) 



And for the parameters of the if-matrix obtain 
IV (892) = 0.0504 GeV, 

Js~ 2 = 1.35 GeV, Js~ z = 1-75 GeV, fa = 0.110, 

fa = -0.685, Ai = 0.142, A 2 = 1.089. 

7^ - 0.204, 7$ = -0.983, 7$ = 8.329, 



tS 5 = 



-0.052, 



7$ = 0.426, 



7$ = -3-834. 

(B8) 

W's in units 



with the 7^°^'s in units of GeV -2 and 7 
of GeV -4 . The phase, (fin and magnitude, |tn| of 
Kit — > Ki t scattering amplitude is compared to the date 



in Fig. 11 We can express i a p = t a p / '(4q a Qp) in terms of 
a product of poles, zeros and the Omnes-Muskhelishvili 
function 



l a/3 



x exp 



N, 



n 



ap- 



is - s (a/3) ) 

Ul=l,s( S ~ S P,0 

ds , <j>ap{s') 
1"i(m K +m w ) 2 s'(s'-s) 



(B9) 



where N z<a p is the number of zeros of t a p for which 
we find N Z}U = ^,22 = 7, N zA2 = N z , 2 i = 6. The 
normalization factors are given by Nn — 7.075, Ni 2 — 
-421.989, N 22 = 2.808. The positions of the poles and 
zeros are given by (in units of GeV 2 ) 

s P1/2 = 0.3573 ± i0.4055, s^ 3/4 = 2.4912 ± i0.5762, 
sp,5/6 = 20.3504 ±i3. 3856, 

s P7 = -0.00489, sp^s = -6.2666, (B10) 
and 

s^] n = 0.3428 ± i0.4470, s?l] A = 2.2188 ± i0.6175, 



2,1/2 
.(") 

'z,5/6 

.(12) _ 
z.1/2 

.(12) 



2,3/4 



12.8061 ± i0.2470, s 



(ii) 
2,7 



= 1.9181 ±z0.3669, s y z 



(12) _ 



3.3956, 



= M 



K* (892)' 



s (U) -0 



.(22) 
2,1/2 

(22) 
'2,5 = 



2.2704 ± i0.2273, sj^ 4 = 20.2775 ± i3.3479, 



-0.00489, s 



(22) 
2,6 



-6.1874, 



.(22) 
'2,7 



(Bll) 



13 




respectively. As can be seen from Fig. [l2j this if-matrix 
leads to a dramatic, most likely unphysical, drop in the 
phase 0n (dashed line) around 6 GeV 2 and this results in 
both (f>n and Sn vanishing asymptotically. Furthermore 
the resulting t- matrix has complex poles and zeros on the 
physical sheet (see Eq. (BIO) and Eq. (Bll )). The origin 
of these unphysical poles can be illustrated by considering 
a single channel, K\% only, with a single pole and constant 
background term. The resulting i-matrix element is then 
given by, 




0.8 1.2 1.6 2 

(GeV) 



FIG. 11: 0ii (top) and [tn| (bottom) of Kn — > Kir scatter- 
ing amplitude vs data from [20] (squares), |21] (circles), and 
[T9] (triangles). 



a 



8(M 2 



if* (892) 



-,)- 



-(m K - 



lT7?l(A4.( 8 92)- 3 ) 

™7r) 2 ][s-(TOjf+TO^) 2 ] ^ J £ 



(B12) 



if* (892) 



If, for simplicity, we replace ttik by and keep only 
the imaginary of Ii(s), in the limit \s\ — > oo and 7^ — > 
with |7ii||s| ^ a\ one finds 



'11 



7n 



1 - 7^711 s 



(B13) 



±- 



the pole is on the first sheet and 
Even though the K-matrix itself 



In the limit 711 
approaches s — ► 
has unphysical singularities and zeros it still faithfully 
reproduces the phase and magnitude of the amplitude of 
Kir scattering data up to 1.8 GeV. Similarly to the cases 
in 7T7T and KK scattering presented in previous sections, 
we will truncate the K-matrix solution at si ow and match 
it with Regge parameterization at Shigh- 



2. Regge parameterization for tt°K ± — ¥ ■k°K ± 

(S > Shigh) 



Asymptotically we only retain Pomeron (P) in the t- 
channel and the K* trajectory in the u-channel 



pi t =o 



(t,s,u) 



1 + e-*™H*> l + ap(t) bt (t) 
jip a P (t) e (s/s) py ' . 



sin Tra.p(t) 



(B14) 



sanra K .(u) H 2a K .(0) + l v K ' 

(B15) 
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The Pomeron trajectory is given in Eq. (All) with pa- 



rameters in Pomeron parameterization are taken from 
[HIM], except the coupling constant /3 p K = fi P) f { / } = 
[/3™px R ]i = 1.709. The K* trajectory is given by 
a K ,(u) = 0.352 + 0.882u as in Section |A~2] From s-t 
and s-u channel crossing, we obtain 



The complete amplitude for the ttK 
given by, 



■kK amplitude is 



tK 



(s) 



4-Kmatrix I \ „ ^ „ 



Regge 
ttK — >7tK 



(S),S > S h igh 



(B18) 



1 



h=0 



Vq Kk 
1 / 



Iu = i 



(u,t,s). 



(B16) 



The P-wave projection of the Regge amplitude in irK 
ttK scattering is given by 



Regge 
1 



(*) = 



(d cos) cos ( 



F 



i t =o 



where we choose ^/si ow — 2.3GeV, y'Shigh — 2.5GeV for 
the real part of the amplitude and ^/Shigh = 2.7GeV 
for the imaginary part of the amplitude as shown in 
Fig. [l2j From the P-wave projection of the Regge am- 
plitude we find the following asymptotic behavior for the 
phase ^>n (s) and denominator function of the complete 
amplitudes, 



i>n — > arctan[- 



sin7rap(0) tt 



1 + cos 7rap(0) 
1 i 



(B19) 



(u, t, s)] . (B17) The phase is shown in Fig. [l2j 
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FIG. 12: Real (left) and imaginary (right) parts of the 
isovector, P-wave amplitude, t^K ->ttK (s) and phase 4>n (solid 
curves). The dashed curves are the result of the if-matrix pa- 
rameterization. 



